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We investigate the possible corrections to the linear Regge trajectories for the light- 
quark meson sector by matching two-point correlators of quark currents to the Operator 
Product Expansion. We find that the allowed modifications to the linear behavior must 
decrease rapidly with the principal quantum number. After fitting the lightest states in 
each channel and certain low-energy constants the whole spectrum for meson masses and 
residues is obtained in a satisfactory agreement with phenomenology. The perturbative 
corrections to our results are discussed. 
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The obs erve d masses squared of mesons with given quantum numbers form linear 
trajectories depending on the number of radial excitation n. This is a strong 
indication that QCD admits an effective string description, as this type of spectrum 
is characteristic e.g. of the bosonic string. In the bosonic string model the slope of 
all trajectories must be equal since this quantity is proportional to the string tension 
depending on gluedynamics only. However, there exist sizeable deviations from the 
string picture. In the present analysis we examine possible corrections to the linear 
trajectories in the vector (V), axial- vector (A), scalar (S), and pseudoscalar (P) 
channels Our method is based on the consideration of the two-point correlators 
of V,A,S,P quark currents in the large- iVc limit of QCD ^. On the one hand, by virtue 
of confinement they are saturated by an infinite set of narrow meson resonances, 
that is, they can be represented by the sum of related meson poles in Euclidean 
space: 

d4xexp(zgx)(grg(a;)grg(0))pia„ar - 02/ 2! v 



expressing the quark- hadron dualitj0. Here J = S, P,V, A; F = i, 75, jfj.,jp,^5- Fur- 
ther we denote Fs.p = Gs.pms.p. On the other hand, their high-energy asymptotics 
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is provided by the perturbation theory and the Operator Product Expansion (OPE) 
with condensates'^. Matching these two approaches results in the Chiral Symmetry 
Restoration sum rules representing a set of constraints on meson mass parameters'^. 
We performed the further analysis in the chiral limit and firstly in the leading order 
of perturbation theory. 

Let us consider the linear ansatz for the meson mass spectra with a non-linear 
correction 5: 

m^j{n)^mlj + an + 5j{n), J = V,A,S,P. (2) 

The last term in Eq. Q signifies a possible deviation from the string picture in QCD. 
The condition of convergence for the generalized Weinberg sum rules ^imposes the 
equality of slopes and a falloff of non-linear correction. 

The asymptotic freedom leads to the relation between residues and masses: 
F^{n) ^ '^"rfi"^ ■ Our analysis showed that the analytical structure of the OPE 
admits, however, exponentially small deviations from this relation (in contrast to 
tolffl), 

(n) ^ Afj + an + A^e-^^-'', (3) 

FlAn) = « + AYe-'^^-^ , G|p(n) = « (y^^ + A^fe-^-^^ (4) 

with certain constants ^^'"^g^ and Bm.F.G > to be fitted. It is plausible to 
suppose that, for masses, the dynamics under these exponentially small corrections 
is governed mostly by gluons and thereby does not depend on flavor. Thus, we keep 
the exponent Bm the same for all channels. For the same reason we regard Bf.g as 
independent of parity. 

In Tables 1 and 2 we show an example of numerical fits resulting from our 
approach where in the S,P-channels two possibilities are considered: 7r-meson 
belongs to the radial Regge trajectory ('7r-in') and does not ('7r-out'). The in- 
puts general for all tables (if any) are: a = (1120 Me V)^, {qq) = -(240 Me V)^, 
ird'-'^f) = (360MeV)^ U = 103 MeV, = 2-^, a, = 0.3. The units are: 
m(n), F(n), G{n) — MeV; A„ — MeV^; Ap, Aa, Bf.g,^ — MeV". 

Let us discuss the impact of running coupling constant at next-to-leading order 
of perturbation theory. We have in this case the contribution to the imaginary part 
of correlator which is related to the full correlator through the dispersion relation. 
In the large- A^c limit for the vector and axial-vector case this leads to: 

lmn{t) = ^ (l + ^) + 0{al) = ^ £ 2F\n)S {t ~ m\n)) . (5) 

^ ^ n=0 

We perform summation in Eq. Q by applying the Euler-Maclaurin summation 
formula. This provides smoothness of this expression. Finally we obtain the relation: 



dn 8tt^ \ IT J 



(6) 
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Table 1 . An example of parameters for the mass spectra. The existing experimental 
values 2>il are displayed in brackets. 



Case 


Inputs 


Fits and constants 


VA 


my(0) = 770 (769.3 ±0.8), 
mA{0) = 1200 (1230 ±40), 


M = 920, B„ = 0.97, Bp = 0.72, 
A]:^ = -500^, A4 = 770^, 
A^, = 0.0012, A^ = -0.0031, 
Lio = -6.5 ■ 10-3, Am^ = 2.3 


SP 
(TT-in) 


ms(0) = 1000, 
mp(0) =0, 

mp(l) = 1300 (1300 ± 100), 
Bm = 0.97 


M = 840, Bo = 0.42, 
AS = 5502, = -8402, 
Ag = -0.0009, Ag = 0.0004, 
Ls = 1.0 ■ 10-3 


SP 

(TT-OUt) 


ms(0) = 1000, 
mp(0) = 1300 (1300 ± 100), 
mp(l) = 1800(1801 ± 13), 
= 0.97 


M = 1470, Bg = 1.27, 
AS = -10802, A^ = -6902, 
Ag, = 0.0213, Ag = 0.0067, 
Lg = 0.9 ■ 10-3 



Table 2. Mass spectrum and residues for the parameter sets of Table 1. 



n out 





1 


2 


3 


mvin) 


770(775.8 ± 0.5) 


1420 (1465 ± 25) 


1820 (1900?) 


2140 (2149 ± 17) 


Fv{n) 


138 (154±8) 


135 


133 


133 


niAin) 


1200 (1230 ± 40) 


1520 (1640 ± 40) 


1850(1971 ± 15) 


2150 (2270 ± 50) 


FAin) 


116 (123 ±25) 


125 


128 


130 


ms{n) 


1000 (980 ± 10) 


1440 (1507 ± 5) 


1800 (1714 ±5) 


2100 


Gs(n) 


176 


178 


178 


179 


mp in) 


(TT-in) 


1300 (1300 ± 100) 


1760 (1801 ± 13) 


2100 (2070 ± 35) 


Gp(n) 




179 


179 


179 


ms{n) 


1000 (980 ± 10) 


1730 (1714 ±5) 


2120 


2420 


GsH 


243 


199 


185 


181 


mpin) 


1300 (1300 ± 100) 


1800 (1812 ±14) 


2150 (2070 ± 35) 


2430 (2360 ± 30) 


Gp{n) 


201 


186 


181 


180 



Eq. © can be approximated by finite differences. For the first two states tlris reads: 



(7) 



TO^, — m^, — rn^^ Stt 

Substituting experimental values into Eq. Q one arrives at the estimate: 1.5 ± 
0.2 « 1.3 ± 0.6 « 1.4 in a good agreement with the phenomenology. 

In order to reproduce the running coupling behaviour in Eq. JSJ one should 
accept the following ansatz for the residues: 

l + 4Uln-^j U exp. corr. (8) 
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However this ansatz introduces into the OPE the additional terms ~ ''^'j (and 
powers of these terms) which are absent in the standart OPE. Thus ansatz ^ 
must be improved. This problem is under our studies now. 
Let us summarize the results of our analysis: 

1) The convergence of the generalized Weinberg sum rules requires the universality 
of slopes and intercepts for parity conjugated trajectories. 

2) The matching to the OPE cannot be achieved by a simple linear parameteriza- 
tion of the mass spectrum, the linear trajectory ansatz. There must exist deviations 
from the linear trajectory ansatz triggered by chiral symmetry breaking. These de- 
viations must decrease at least exponentially with n. 

3) For heavy states, the D-wave vector mesons have to decouple from asymptotic 
sum rules. This fact implies the exponential (or faster) decreasing the corresponding 
decay constants F|,(n). 

4) Our results seem to exclude a light ct(600) particle as a quarkonium state and 
rather favor the non-linear realization of chiral symmetry with the lightest scalar of 
mass ~ 1 GeV, its chiral partner being the 7r'(1300). 

5) In our approach the quantities LgjLio and AttItt are obtained, in satisfactory 
agreement with the phenomenology. 

6) Perturbative corrections can be systematically treated making, nevertheless, 
small effects on the fits presented. 
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